Curvature induced toroidal bound states 

Mario Encinosa and Lonnie Mott 
Department of Physics 
Florida A & M University 
Tallahassee, Florida 32307 

Abstract 

Curvature induced bound state {E < 0) eigenvalues and eigenfunctions for a particle 
constrained to move on the surface of a torus are calculated. A limit on the number of 
bound states a torus with minor radius a and major radius R can support is obtained. 
' A condition for mapping constrained particle wave functions on the torus into free 

particle wave functions is established. 

Pacs number(s): 03.65Ge, 68.65.-k 
1. Introduction. 
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^ ■ The physics of nanostructures [1,2] and quantum waveguides [3-7] may make questions 

concerning curved surfaces in quantum theory increasingly relevant to device modelling. 
Many workers have investigated with varying levels of formal machinery the existence of 
bound states of quantum systems on curved strips, tubes and hypervolumes [8-12]. The 
00 ■ common thread through much of the work described in [8-12] is the existence of an attractive 

^ ■ potential that appears in the Schrodinger equation as a consequence of constraining a particle 

O ■ from higher to lower dimensional manifolds. (In the majority of work the dimensionality is 

■ reduced from three to two, but see [11] for the generalization to other cases.) This potential, 
called here the curvature potential Vc, has been shown sufficient to cause bound states in 
model systems. In this work Vc for a particle constrained to the surface of a torus is derived 

^ ■ and bound state surface toroidal wavefunctions (STWs) are calculated. 

^ . This brief report is organized as follows: in section 2 the method by which Vc is derived 

O^' is concisely described for a symmetric but non-trivial geometry. The method is then applied 

to the torus. In section 3 a brief description of the procedure used to solve the Hamiltonian 
^ . found in section 2 is given and some low-lying bound state eigenvalues and STWs shown. 

■ Conclusions appear in section 4. 

2. Derivation of Vq. 

In the interest of clarity the derivation of Vc will be performed for a cylindrically sym- 
metric surface. The extension to the general case is straightforward and the salient points 
still obtain. 

Let ep, e^, ez be standard cylindrical coordinate system unit vectors. A cylindrically 
symmetric surface may be described by the Monge form 

r{p,(f)) = p ep + S{p) e,. (1) 

S{p) gives the shape of the surface. Points near the surface S{p) may be described by 

x(p, 0, q) = r(p, 0) + g e„ (2) 
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with e„ everywhere normal to the surface. The metric near the surface S is 
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d(l)'^+ dq' 



= + qk^Y dp'' + + qk^f d(l>' + dq' 



(3) 
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with subscripts indicating differentiation and Z 
straightforwardly from 



1 + 5"^. The Laplacian can be found 
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dq^ 
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but there is no advantage to writing the the Laplacian explicitly until the constraint that 
places the particle on the surface is effected. 

Consider a situation where a large confining potential everywhere normal to S acts to 
restrict the particle to S. This potential, called Vn{q) here, could take a hard wall or oscillator 
form, but however chosen it causes — > 0. In this limit the wave function is expected to 
decouple into surface and normal parts, or in the language of [13], into "fast" and "slow" 
functions 

^ Xs(P»Xn(g)- (6) 

Conservation of the norm in the decoupled limit implies [11,14,15] 



l^afWdSdq = \xs\^\xnfdSdq 



(7) 



where 1^ = 1 + 2qH + q'K and dS the surface measure. H, K are the mean and Gaussian 
curvatures given by 

(8) 



Write 



H^-{k, + k2), 

K — kik2- 

^ _ XsXn 



(9) 
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and insert the right hand side of equation (10) into the time independent Schrodinger equa- 
tion. Performing the differentiations and taking g — > gives the pair of equations 

Id' Id Z.d Id' 



Z^dp'^ Zpdp Z^dp p'^dcj)'^ 
Id'xn 



Xs 



EsXs 



2 dq 



,2 + K(g)Xn = EnXn- 



The curvature potential Vc is 



Vc^--[H'-K\. 



(11) 
(12) 

(13) 



The surface given by equation (1) was chosen to illustrate that there are modifications to 
the Laplacian (aside from the appearance of Vc) even for surfaces possessing symmetry. In 
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[16] it was shown that kinetic energy modifications for some parameterizations of S{p) can 
be small, but it is easy to conceive of cases where the contrary would be true. 

Now apply the above procedure to the torus. Let F — R+a cos9. Points near the surface 
of the torus may be parameterized as 

:x.{6,(j),q) = Fep + a sinO ez + q Bn- (14) 
Proceeding as above gives the constraint and curvature modified Hamiltonian 

(15) 



1 sine d 1 W 1 

-Ho-^^,^Ho + Vc. (16) 
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In equation (16), Hq is identical to — derived from eq. (14) with g = 0; it is the 
Hamiltonian for particle on a toroidal surface subject to no other potential [17] (it proves 
convenient to refer to the zero potential case as the free system) . In contrast to the operator 
that appears in equation (11), no surface dependent prefactors, i.e. Z{p) or additional terms 
modifying Hq are present. It is interesting to compare equation (15) to the corresponding 
operator for spherical and cylindrical surfaces; for those surfaces the constraint proccdiire 
gives a coordinate independent Vc behaving as ~ —l/B? [18] and again no modifications to 
Hq. This point will be discussed further in section 4. 

3. Solution method; results. 

Setting a — (5 — 2Ea'^ and making the standard ansatz for the azimuthal eigenfunction 
x(0) = exp [imcf)] in equation (15) gives 



9^-0 a sine dijj (m^a^-i) 



+ 0. (17) 



[i + a cos e] de [1 + a cos e\ 

Equation (17) can be solved numerically, but it is convenient to have approximate analytic 
representations of its eigenfunctions. Recently a method was found for obtaining closed form 
solutions for zero energy states and for surface potentials Vs{0) which satisfy an auxiliary 
condition derivable from Hq [19]. However, the Vs{0) — case is no more easily solved with 
the method given in [19] than the method employed below. 

Solutions of equation (17) can be found by defining z — exp[i9] and writing 

oo 

i^iz)^ ^ c„^^ (18) 

n=—oo 

The Hamiltonian given by equation (15) is invariant under 9 — > —9, so the solutions of 
equation (17) can be split into odd and even parity eigenfunctions, yielding a scries in sines 
for negative parity states and cosines for positive parity states. Computing the eigenvalues 
and eigenfunctions of equation (17) follows from a method given in detail in [17]. 
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Tabic 1 shows eigenvalues and wave functions for a = 3/4,1/2,1/4,1/20 for those m 
values which yield (3 < states (only three states are given for a = 1/20; there are nine 
total). No negative parity states appear in table 1. In table 2 wave functions for states 
corresponding to those in table 1 with Vc shut off are given. As evidenced in table 2, the 
ground state m — wave function of Hq for any a is a constant. A natural question is: 
Should the Vc ^ ground state wave function be compared to the constant wave function or 
the lowest /5 7^ wave function? Here the constant wave function was chosen on the grounds 
that it is the state actually altered from its constant value by Vc- 

4. Conclusions. 

In this brief report wave functions for bound states of a particle constrained to the surface 
of a torus were obtained for several values oi a — a/ R. Constraint and curvature effects were 
shown to alter the angular dependence of the free particle STWs. 

An interesting consequence of equation (17) is manifested by the results presented in 
table 1. For m ^ 0, a = ^ provides a cutoff for the existence of bound states. It follows 
that there are no m 7^ bound states for a > 1/2. Additionally a = l/2m provides a series 
of magic radii for which an m 7^ state of the constrained system maps exactly into the 
m = state with the same n and parity of the Vc = free STW. 

In section 2 it was stated that the torus shares with the sphere the property that constraint 
adds only a curvature potential to the Hamiltonian, leaving on the surface unchanged. 
This comparable behavior is likely a consequence of the torus being the most symmetric 
compact genus one surface that can be embedded in R^. It would be interesting to learn if 
there are higher genus surfaces embedded in R^ for which H — — reduces to the lower 
dimensional operator Hq plus a curvature potential upon imposing the condition given by 
equation (10). This question generalizes to M < N dimensional surfaces embedded in R^ . 
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Table 1: Eigenvalues and wave functions of equation (17) for four values of a. Normalizations 
(not including the (27r)~2 from the (j) dependence) are in brackets proceeding the functions. 
Terms not shown are at least an order of magnitude smaller than those hsted. There are no 
negative parity states with /3 < 0. 



a 


/5 




(^); VcT^o 






.75 


-1.0725 


*10 


= (. 1298) [4.6072- 


5.2143 cos^ + 2.2465 cos2^ - .9495 cosW] 


.50 


-0.3512 


*10 


= (.2455) [2.4509- 


.9015 cos^ + 


.1921 cos2^] 


.25 


-0.2673 


^10 


= (.3765) [2. 1458- 


.2916cos^ + 


.0280cos2^] 


.25 


-0.1987 


^11 


= (.3826)[2.1069- 


.2138cos^ + 


.0197cos2^] 


.05 


-0.2506 


^10 


= (.8813) [2.0254- 


.0508cos^] 




.05 


-0.2481 


*11 


= (.8814)[2.0251 - 


.0507cos^] 




.05 


-0.2406 


*12 


= (.8817) [2.0244- 


.0487cos^] 





Table 2: Eigenvalues and wave functions of Hq corresponding to those appearing in table 1. 



a 


/3 


^ nn 


m vc^Q 




.75 


0.0000 


^00 


= .4607 




.50 


0.0000 


^00 


= .5642 




.25 


0.0000 


^00 


= .7979 




.25 


0.0641 


^'ll 


= (.4073) [1.9676 + 


.0648cose] 


.05 


0.0000 




= 1.7841 




.05 


0.0025 


^11 


= (.8822) [1.9998 + 


.0005cos^] 


.05 


0.0010 


*12 


= (.8822) [1.9996 + 


.0002cos^] 
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